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Executive  Summary 

Introduction 

Software  has  been  written  in  order  to  provide  U.S.  Army  Research 
Laboratory  (ARL)  scientists  prototyping  a  Meteorological  Measuring  Set 
-  Profiler  (Proof-of-Concept)  (MMS-P  (POC))  System  at  White  Sands 
Missile  Range  with  a  "real-time"  on-site  tool  to  investigate  the  quality  of, 
in  particular,  meteorological  messages  output  by  the  MMS-P  (POC) 
mesoscale  model.  The  software  consists  of  Bourne  shell  scripts, 
FORTRAN,  C,  and  GNUPLOTt  code.  There  are  seven  shell  scripts  that 
call  on  the  FORTRAN  and  C  code  modules  for  data  reformatting  and 
processing.  The  GNUPLOT  code  modules  allow  for  either  the  visual 
display  of  results  or  the  creation  of  a  hardcopy  of  the  graphical  results. 
The  graphical  results  contain  the  difference  between  the  truth  data  and 
the  data  being  analyzed  for  these  variables:  pressure,  virtual  temperature, 
u-  and  v-component  of  the  wind. 

Overview 

Six  different  types  of  comparisons  using  varying  formats  of 
meteorological  data  can  be  made.  For  the  MMS-P  (POC)  testing,  "truth" 
data  is  derived  from  two  sources.  They  are  (1)  the  Meteorological 
Measuring  Set  (MMS),  which  is  an  Army  truck  instrumented  to  collect 
meteorological  data  from  a  balloon  launch.  The  data  collected  from  these 
launches  can  be  in  either  computer  meteorological  message  (metcm) 
format  or  in  a  format  similar  to  those  from  the  National  Weather  Service 
(NWS).  (2)  NWS  radiosonde  observation  (raob)  collected  from  the 
Internet.  The  data  scrutinized  for  accuracy  during  the  MMS-P  (POC)  test 
were  predominately  the  metcm's  output  by  the  Computer-Assisted 
Artillery  Meteorology  Battlescale  Forecast  Model  (CAAM  BFM) 
mesoscale  model  but  occasionally  the  CAAM  BFM  binary  output  file 
containing  forecasted  variables  was  compared.  Other  types  of 
comparisons  can  be  made  as  well.  The  user  has  the  options  listed  below. 

1.  MMS  balloon  data  versus  CAAM  BFM  forecast  data  (binary  file) 

2.  MMS  balloon  data  versus  NWS  format  balloon  data 

3.  NWS  format  balloon  data  versus  CAAM  BFM  forecast  data 


t  Freeware  graphics  display  package 
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4.  NWS  format  balloon  data  versus  like  data  (temporally  and/or 
spatially  displaced) 

5.  CAAM  BFM  (computer  meteorological  message  [metcm]);  American 
Standard  Code  for  Information  Interchange  [ASCII]  file)  versus  MMS 
metcm 

6.  CAAM  BFM  metcm  versus  like  data  (temporally  and/or  spatially 
displaced) 

All  of  the  interpolation,  both  horizontal  and  vertical  takes  place  within 
the  FORTRAN  modules.  The  C  modules  handle  the  parsing  of  data  files 
and  the  derivation  of  statistics. 

Conclusions 

Based  on  preliminary  testing,  this  software  has  proven  to  be  beneficial  in 
determining  the  quality  of  meteorological  data  output  by  CAAM  BFM 
both  during  and  after  MMS-P  (POC)  testing. 

A  new  method  for  the  display  and  statistical  analysis  of  several 
meteorological  data  types  has  been  developed  by  the  creation  of  new 
software  in  both  FORTRAN  and  C  and  the  use  of  existing  FORTRAN 
modules  as  well  as  the  GNUPLOT  freeware  display  program.  The 
impetus  for  this  software  development  was  the  need  for  real-time  on-site 
statistical  analysis  at  a  follow-on  test  of  the  MMS-P  (POC)  by  ARL 
scientists  at  White  Sands  Missile  Range,  NM  9-20  August  1999.  The 
MMS  as  it  is  currently  fielded  provides  meteorological  data  for  artillery 
battery  aiming  adjustments.  Its  data  is  derived  from  a  weather  balloon 
release  and  the  subsequent  tracking  of  the  balloon  and  recording  of  the 
data.  The  next  generation  MMS,  the  MMS-P,  is  being  developed  so  that 
rapid  updates  of  meteorological  information  can  occur  with  obvious 
benefit  to  artillery  batteries  that  will  receive  the  most  current  data  for 
aiming  adjustments.  In  addition,  it  is  hoped  that  eventually  the  weather 
balloon  can  be  removed  from  the  battlefield  and  can  be  replaced  by  a 
profiling  suite  of  equipment. 


1.0  Overview 

There  were  two  key  parts  to  this  Meteorological  Measuring  Set  -  Profiler 
(Proof-of-Concept)  MMS-P  (POC)  test.  The  first  is  to  verify  the 
interoperability  of  all  hardware  components  including  a  radar, 

.  radiometer,  surface  sensor,  satellite  receiver,  as  well  as  several  computer 
servers  for  data  communication.  The  second  is  to  ascertain  the  accuracy 
of  the  meteorological  mesoscale  forecasting  model  output  when  profiler 
data  from  MMS-P  (POC)  is  used  as  one  of  the  inputs.  The  model  for 
MMS-P  (POC)  is  the  Computer-Assisted  Artillery  Meteorology 
Battlescale  Forecast  Model  (CAAM  BFM).  It  is  a  hydrostatic  model 
capable  of  providing  12-hour  forecasts  and  yields  two  outputs  that  are  of 
relevance  to  this  software  development.  They  are  the  computer 
meteorological  message  (metcm)  and  the  binary  output  (out.bin)  file, 
which  contains,  among  other  things,  the  u-  and  v-components  of  the 
wind,  ambient  temperature,  and  mixing  ratio.  The  metcm  is  specifically 
tailored  for  the  needs  of  the  artillery  and  is  a  27-line  message  containing 
data  from  the  surface  up  to  20-km  above  ground  level  (AGL).  At  each 
metcm  level  are  the  level  number  and  four  variables:  wind  direction  (tens 
of  mils),  wind  speed  (kt),  virtual  temperature  (K),  and  pressure  (mb).  The 
mil  is  a  unit  particular  to  artillery.  To  convert  from  mils  to  degrees  one 
multiplies  by  the  factor,  360°/ 6400  mils. 

In  order  to  determine  the  accuracy  of  CAAM  BFM  output,  some  data  sets 
had  to  be  chosen  as  "truth".  For  this  particular  MMS-P  (POC)  test  at 
White  Sands  Missile  Range,  Meteorological  Measuring  Set  (MMS) 
weather  balloons  were  launched  quite  often  in  the  proximity  of  the 
MMS-P  (POC)  equipment.  Also,  data  from  the  weather  balloon  launched 
at  Santa  Teresa,  NM  (near  El  Paso,  TX)  was  available  for  all  but  one  day 
of  the  experiment  soon  after  launches  at  0000  Universal  Time  Coordinates 
(UTC)  and  1200  UTC.  Also  noteworthy  is  that  within  each  MMS  message 
the  data  can  be  found  in  two  different  formats  and  are  not  identical  data 
sets.  One  format  is  the  same  type  used  by  the  National  Weather  Service 
(NWS)  that  reports  "point"  measurements  of  temperature,  wind,  and 
moisture  parameters  at  the  so-called  "mandatory"  levels  plus 
"significant"  levels.  The  other  is  the  metcm  where  the  wind  and 
temperature  are  reported  as  an  average  over  a  particular  height  range 
with  pressure  being  the  only  reported  point  measurement.  In  addition, 
the  format  of  the  NWS  weather  balloon  data  posted  at  various  Internet 
sites  varies  somewhat.  Specifically,  sometimes  the  reports  contain  a  line 
listing  the  station  elevation  and  sometimes  this  is  absent.  In  addition. 
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some  reports  will  occasionally  report  pressure  to  a  tenth  of  a  millibar 
rather  than  as  a  simple  integer.  Thus,  the  analysis  software  was  designed 
to  handle  these  variable  formats  and  to  be  flexible  enough  to  allow  the 
user  to  compare  any  of  the  data  types  discussed  above.  In  addition,  so 
that  the  user  could  easily  visualize  the  results  of  a  given  comparison, 

GNUPLOT$  was  adopted,  primarily  due  to  the  author's  familiarity 
with  it. 
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2.0  Statistical  Analysis 

Because  the  CAAM  BFM  data  is  gridded  data,  bilinear  interpolation  is  a 
practical  means  of  interpolating  horizontally  to  a  particular  point,  in  this 
case,  to  the  latitude/longitude  of  the  MMS-P  (POC)  experimental  site  at 
White  Sands  Missile  Range.  For  8  of  the  12  options  available  to  the  user, 
vertical  interpolation  to  fixed  heights  mean  sea  level  (MSL)  is  carried  out 
so  those  variables  from  two  vertical  profiles  are  compared  at  the  same 
heights.  There  are  only  six  distinct  types  actually;  however,  since  there 
are  print/display  options  for  each  type,  the  user  will  have  12  choices. 
The  four  types  which  involve  vertical  interpolation  are  comparisons 
between:  1)  MMS  radiosonde  observation  (raob),  model  2)  MMS  raob, 
NWS  format  raob  3)  NWS  format  raob,  model  and  4)  2  NWS  format 
raobs.  Simple  linear  interpolation  in  the  vertical  can  be  carried  out  for  the 
u-  and  v-components  and  for  virtual  temperature.  However,  for 
pressure,  the  hypsometric  relationship  (which  states  essentially  that 
pressure  decreases  exponentially  with  height)  is  used.  This  horizontal 
and  vertical  interpolation  is  carried  out  in  three  FORTRAN  programs. 

1.  compare_raobl_raob2.f 

For  the  case  of  a  comparison  between  two  NWS  format  raobs. 

2.  compare_raob_BFMout_bin.f 

For  the  case  of  NWS  format  raob  and  model  data  comparison. 

3.  compare_raob_MMSraob.f 

For  the  case  of  MMS  raob  and  NWS  format  raob  comparison. 

When  a  MMS  raob  is  used,  the  C  program  parseMMS.c  must  be  called. 
MMS  messages  contain  both  a  vertical  profile  like  that  from  a  raob  from 
the  NWS  and  the  same  data  in  a  metcm  format  plus  lots  of  other 
extraneous  information.  parseMMS  looks  for  and  extracts  the  NWS 
format-sounding  portion  of  the  MMS  message.  Similarly,  when  a  raob 
such  as  those  from  the  NWS  is  used,  the  C  program,  parselntemetRaob.c 
extracts  the  appropriate  data. 

The  C  program  rd_raob_compare.c,  using  the  output  from  these 
FORTRAN  programs,  then  checks  for  missing  values  and  calculates  the 
root  mean  square  error  (rmse)  for  the  u-  and  v-components,  virtual 
temperature,  and  pressure. 
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There  are  two  comparison  types  in  which  no  vertical  interpolation  is 
done: 

1.  CAAM  BFM  metcm,  MMS  metcm 

2.  CAAM  BFM  metcm,  CAAM  BFM  metcm. 

For  CAAM  BFM  metcm's,  the  C  program, 
reformat_bfm_metcm_to_columns  is  called.  For  MMS  metcm's,  the  C 
program,  reformat_mms_metcm_to_columns  is  called.  The  latter  is  much 
more  complex  as  it  must  search  for  patterns  within  the  MMS  message  to 
extract  only  the  metcm  portion.  The  deltas  between  variables  and  rmse 
values  for  u,  v.  Tv  and  P  are  then  calculated  by  the  C  programs, 
compare_mms_bfm_metcm  or  compare_bfml_bfm2_metcm  depending 
on  the  option  selection.  It  is  emphasized  that  if  one  is  going  to  use  either 
of  these  two  options  that  the  two  stations  being  compared  have  only 
minor  surface  elevation  differences  and  preferably  no  major  topographic 
features  intervening  such  as  mountains. 

Seven  scripts  control  the  entire  statistical  analysis  and  display /print 
process.  The  primary  script  is  called  "stats"  and  it  in  turn  calls  one  of  the 
other  six  scripts  depending  on  user  response  to  queries.  The  six  types  of 
comparisons  that  the  user  has  to  choose  from  and  the  respective 
controlling  scripts  are  listed  below. 

1.  MMS  balloon  data  compared  to  CAAM  BFM  out.bin  data 
MMSraob_model.sh 

2.  MMS  balloon  data  compared  to  NWS  format  balloon  data 
MMSraob_raob.sh 

3.  NWS  format  balloon  data  compared  to  CAAM  BFM  out.bin  data 
raob_model.sh 

4.  NWS  format  balloon  data  compared  to  NWS  format  balloon  data  for 
another  site/time 

raob_raob.sh 

5.  A  CAAM  BFM  metcm  compared  to  a  MMS  metcm 
caam-metcm_mms-metcm.sh 

6.  A  CAAM  BFM  metcm  for  one  location  compared  to  a  CAAM  BFM 
metcm  for  another  location 

bfm-metcm_bfm-metcm.sh 

For  each  of  these  six  cases,  the  user  has  the  option  of  viewing  the 
graphical  results  through  GNUPLOT  or  sending  the  graphical  results  to 


the  computer's  default  printer  (which  must  have  postscript  printing 
capability)  after  a  postscript  file  has  been  created  (by  GNUPLOT).  When 
the  user  chooses  the  display  option,  they  will  see  popping  up  on  their 
screen  respectively,  plots  of  the  u-component  of  the  wind,  v-component 
of  the  wind,  virtual  temperature,  and  pressure.  Each  plot  will  display  for 
15  seconds  before  giving  way  to  the  next  plot.  The  time  each  display  is  on 
the  screen  is  controlled  in  C  code  that  makes  GNUPLOT  calls  (see 
Appendix  C  for  an  example).  The  GNUPLOT  files  are  set  to  "autoscale" 
in  both  the  x-  and  y-directions.  When  the  user  chooses  the  print  option,  a 
postscript  file  is  created  and  automatically  printed  to  the  default  printer 
for  the  computer  they  are  utilizing  (thus  the  user  must  be  cognizant  of 
which  printer  that  is).  In  addition,  regardless  of  the  option  chosen,  the 
root  mean  square  error  for  the  u-  and  v-components  of  the  wind,  virtual 
temperature,  and  pressure  will  be  printed  to  the  screen. 

Figure  1  is  the  initial  stats  menu. 


Figure  1.  Initial 
stats  menu. 


Make  a  selection  from  these  choices: 

1 .  Compare  MMS  raob  uith  BFM  out. bin  data  and  display  results 

2.  Compare  MMS  raob  with  BFM  out. bin  data  and  send  results  to  printer 

3.  Compare  MMS  raob  with  a  raob  and  display  results 

A.  Compare  MMS  raob  uith  a  raob  and  send  results  to  printer 

5.  Compare  a  raob  uith  BFM  out. bin  data  and  display  results 

6.  Compare  a  raob  uith  BFM  out. bin  data  and  send  results  to  printer 

7.  Compare  a  raob  uith  another  raob  and  display  results 

8.  Compare  a  raob  uith  another  raob  and  send  results  to  printer 

9.  Compare  a  CAAMBFM  metcm  to  a  MMS  metcm  and  display  results 

10.  Compare  a  CAAMBFM  metcm  to  a  MMS  metcm  and  send  results  to  printer 

11.  Compare  2  CAAMBFM  metcm 's  and  display  results 

12.  Compare  2  CAAMBFM  metcm's  and  send  results  to  printer 

13.  Quit  the  analysis  and  display /print  program 
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Error  handling  is  built  into  the  stats  script.  For  example,  if  the  user 
selects  outside  the  range  1-13  they  will  see  a  message  indicating  the 
problem  and  how  to  rectify  it.  The  following  graphic  (figure  2)  is  an 
example  of  an  invalid  choice. 


Figure  2.  Result  of 
user  making  an 
invalid  selection  in 
stats. 


10.  Compare  a.  CAAMBFM  met  cm  to  a.  MMS  met  cm  and  send  results  to  printer 
11-  Compare  2  CAAMBFM  met cm's  and  display  results 

12.  Compare  2  CAAMBFM  met cm’s  and  send  results  to  printer 

13.  Quit  the  analysis  and  display/print  program 


6  is  not  a  valid  choice  (must  he  in  range  1-13) 

Make  a  selection  from  these  choices: 

1-  Compare  MMS  raob  with  BFM  out  .bin  data  and  display  results 

2.  Compare  MMS  raob  with  BFM  out  .bin  data  and  send  results  to  printer 

3.  Compare  MMS  raob  with  a  raob  and  display  results 

4.  Compare  MMS  raob  with  a  raob  and  send  results  to  printer 

5.  Compare  a  raob  with  BFM  out .bin  data  and  display  results 

b.  Compare  a  raob  with  BFM  out .bin  data  and  send  results  to  printer 

7.  Compare  a  raob  with  another  raob  and  display  results 

8.  Compare  a  raob  with  another  raob  and  send  results  to  printer 

9.  Compare  a  CAAMBFM  met cm  to  a  MMS  met cm  and  display  results 

10.  Compare  a  CAAMBFM  met cm  to  a  MMS  met cm  and  send  results  to  printer 

11.  Compare  2  CAAMBFM  met cm’s  and  display  results 

12.  Compare  2  CAAMBFM  met cm’s  and  send  results  to  printer 

13.  Quit  the  analysis  and  di splay /print  program 


Figure  3.  Option  1 
query. 


Figure  3  shows  what  the  user  will  see  when  they  select  option  1  (or 
option  2).  Since  the  CAAM  BFM  mesoscale  model  always  outputs  the  file 
generically  named  out.bin,  this  is  the  name  the  stats  program  expects  to 
see  when  model  comparisons  are  to  be  done.  What  the  user  will  need  to 
know  about  the  model  run  is  how  many  hours  of  forecast  were  done. 
They  then  will  be  able  to  index  to  the  particular  forecast  hour  of  interest 
by  giving  that  value  as  in  figure  3.  Here  the  user  has  chosen  the  first  hour 
of  the  forecast.  In  addition,  the  user  will  be  asked  for  a  MMS  balloon  file 
name.  (Note:  the  reader  will  see  in  upcoming  graphics  the  phrase  "raob", 
in  other  words  the  balloon  data  is  picked  up  at  the  ground  via  telemetry.) 


Make  a  selection  from  these  choices: 

1.  Compare  MMS  raob  with  BFM  out.bin  data  and  display  results 

2.  Compare  MMS  raob  with  BFM  out.bin  data  and  send  results  to  printer 

3.  Compare  MMS  raob  with  a  raob  and  display  results 

4.  Compare  MMS  raob  with  a  raob  and  send  results  to  printer 

5.  Compare  a  raob  with  BFM  out.bin  data  and  display  results 

6.  Compare  a  raob  with  BFM  out.bin  data  and  send  results  to  printer 

7 .  Compare  a  raob  with  another  raob  and  display  results  \ 

8.  Compare  a  raob  with  another  raob  and  send  results  to  printer 

9.  Compare  a  CAAMBFM  metcm  to  a  MMS  metcm  and  display  results 

10.  Compare  a  CAAMBFM  metcm  to  a  MMS  metcm  and  send  results  to  printer 

11.  Compare  2  CAAMBFM  metcm^s  and  display  results 

12.  Compare  2  CAAMBFM  metcm "s  and  send  results  to  printer 

13.  Quit  the  analysis  and  display /print  program 

1 

enter  MMS  raob  file  name  (the  BFM  file  is  assumed  to  be  named  out.bin) 
99081616.txt 

enter  the  forecast  hour  you  want 
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3.0  Case  Results 


3.1  Comparison  of  MMS  Balloon  Data  with  CAAM  BFM 
Model  Data  at  a  Particular  Latitude /Longitude 

For  figures  4-7,  the  MMS  data  is  from  a  balloon  launch  at  the  MMS-P 
(POC)  profiler  test  site  on  16  August  1999  at  2200  UTC.  The  model  run  is 
valid  for  16  August  1999, 2000  UTC.  In  the  u-component  comparison,  one 
can  see  that  in  the  lower  2  km,  that  the  model  indicates  a  weaker  westerly 
component.  At  5  km,  the  situation  is  reversed.  For  the  v-component,  the 
graph  indicates  the  model  has  a  stronger  southerly  component  near  the 
surface,  by  about  3  m/s  with  the  situation  reversing  again  at  higher 
altitudes.  For  the  virtual  temperature  comparison  in  figure  6,  the  model 
is  overestimating  surface  temperature  by  almost  4°  K.  This 
overestimation  of  temperature  near  the  surface  is  a  known  current 
CAAM  BFM  model  problem  and  is  being  investigated.  For  pressure 
there  is  fairly  good  agreement  at  all  levels  except  at  1  km  and  above  7  km. 
In  the  BFM  model,  the  levels  for  pressure  and  temperature  are  distinct. 
For  each  pressure  interpolation  level,  the  pressure  at  the  level  below  and 
the  mean  temperature  for  the  layer  are  involved.  As  these  two  variables 
are  defined  at  different  height  levels,  the  interpolation  becomes  quite 
convoluted.  It  may  be  possible  to  improve  the  pressure  interpolation 
slightly.  Dr.  Henmi,  ARL,  has  suggested  considering  interpolating  the 
pressure  downward  and/or  taking  the  pressure  change  with  height  at  the 
center  of  each  grid  point  as  the  average  of  the  pressure  change  with 
height  at  the  four  surrounding  grid  points. 
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Figure  4.  MMS  raob 
for  Profiler  site  (16 
August  1999,  2200 
UTC),  model 
"profile"  for 
Profiler  site  (valid 
16  August  1999, 
2000  UTC)  u- 
component 
deltas  (m/s). 


Figure  5.  MMS  raob 
for  Profiler  site  (16 
August  1999,  2200 
UTC),  model 
"profile"  for  Profiler 
site  (valid  16  August 
1999,  2000  UTC)  v- 
component 
deltas  (m/s). 


Figure  6.  MMS  raob 
for  Profiler  site  (16 
August  1999,  2200 
UTC),  model 
"profile"  for 
Profiler  site  (valid 
16  August  1999, 
2000  UTC)  Tv  deltas 
(deg  K). 


Figure  7.  MMS  raob 
at  Profiler  site  (16 
August  1999, 2200 
UTC),  model 
"profile"  for  Profiler 
site  (valid  16  August 
1999,  2000  UTC)  P 
deltas  (mb). 
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3.2  Comparison  of  MMS  Balloon  Data  with  NWS  Format 
Balloon  Data 

In  this  case  vertical  interpolation  to  fixed  heights  mean  sea  level  (MSL)  is 
carried  out  for  all  the  variables  (u,  v.  Tv,  and  P)  for  both  the  MMS  balloon 
and  the  NWS  format  balloon  data. 


Figure  8  shows  the  u-component  differences  for  the  NWS  raob  and 
MMS  raob.  At  the  surface  the  Santa  Teresa  raob,  which  was  taken 
2  h  later  than  the  MMS  raob,  shows  stronger  westerlies  at  the 
surface. 


Figure  8.  MMS 
raob  at  Profiler  site 
(10  August  1999, 
2200  UTC)  and 
Santa  Teresa,  NM 
raob  (11  August 
1999, 0000  UTC)  u- 
component 
comparison  (m/ s). 


18 


Figure  9.  MMS 
raob  at  Profiler  site 
(10  August  1999, 
2200  UTC)  and 
Santa  Teresa,  NM 
raob  (11  August 
1999,  0000  UTC)  v- 
component 
comparison  (m/ s). 


The  v-component  differences  are  shown  in  Figure  9.  The  MMS  raob 
indicates  a  stronger  southerly  component  at  the  surface.  Above  1.5  km, 
the  two  systems  agree  fairly  well. 
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Figure  10  shows  the  virtual  temperature  differences.  The  earlier  launch 
of  the  MMS  raob  (1600  local  time)  accounts  for  the  warmer  MMS  surface 
virtual  temperature.  Overall,  the  differences  are  small. 


Figure  10.  MMS 
raob  at  Profiler  site 
(10  August  1999, 
2200  UTC)  and 
Santa  Teresa,  NM 
raob  (11  August 
1999, 0000  UTC)  Tv 
comparison  (deg  K). 
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Pressure  differences  are  displayed  in  figure  11.  At  the  surface,  the  raob 
pressure  is  slightly  higher  which  makes  sense  as  part  of  the  diurnal 
pressure  cycle.  As  the  temperatures  cool,  air  density  increases  leading  to 
an  increase  in  pressure.  A  fraction  of  the  difference  is  also  attributable  to 
the  raob  launch  site  being  about  13  m  lower  than  the  MMS  launch  site. 


Figure  11.  MMS  raob 
at  Profiler  site  (10 
August  1999,  2200 
UTC)  and  Santa 
Teresa,  NM  raob  (11 
August  1999,  0000 
UTC)  P  comparison 
(mb). 
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3.3  Comparisons  of  NWS  Format  Balloon  Data  with  like 
Data  Displaced  Either  Spatially  or  Temporally  or  Both. 

Figures  12-15  show  the  deltas  for  u,  v,  Tv,  and  P  for  9  and 
10  August  1999,  1200  UTC  raobs  taken  at  Santa  Teresa,  NM.  Note  in 
figure  13,  the  10  August  1999  raob  indicates  significantly  stronger 
southerly  flow.  Above  2  km,  the  two  raobs  show  only  minor  differences 
as  one  typically  gets  under  the  weak  flow  of  a  ridge  in  summer.  Figure 
14  shows  the  10  August  raob  indicating  a  temperature  2.5  C  °  K  warmer 
at  the  surface  and  2°  K  cooler  at  2  km.  Above  that,  agreement  is 
good — roughly  within  a  degree. 
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Figure  15.  P  deltas  (mb)  for 
two  1200  UTC  Santa  Teresa, 
NM  raobs,  9  August  and  10 
August  1999. 


3.4  Comparison  of  a  CAAM  BFM  metcm  (from  a  Model 
Rim)  with  a  MMS  metcm 

Figures  16-19  are  graphical  comparisons  of  the  differences  between  a 
CAAM  BFM  model  computation  of  u,  v,  Tv,  and  P  as  output  in  a  metcm 
valid  at  Profiler  site  and  the  measurement  of  u,  v.  Tv,  and  P  by  a  MMS 
balloon  output  in  metcm  format  also  valid  at  Profiler  site. 


Figure  16.  U-component 
deltas  (m/s)  between 
CAAM  metcm  for  Profiler 
site  (valid  10  August  1999, 
1200  UTC  and  MMS  metcm 
(10  August  1999,  2200  UTC) 
for  Profiler  site. 
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Figure  17.  V- 
component  deltas 
(m/ s)  between 
CAAM  metcm  for 
Profiler  site  (10 
August  1999, 1200 
UTC)  and  MMS 
metcm  (10  August 
1999,  2200  UTC)  for 
Profiler  site  . 


Figure  18.  Tv  deltas 
(deg  K)  between 
CAAM  metcm  for 
Profiler  site  (10 
August  1999, 1200 
UTC)  and  MMS 
metcm  (10  August 
1999,  2200  UTC)  for 
Profiler  site. 


Figure  19.  P  deltas 
(mb)  between 
CAAM  metcm  for 
Profiler  site  (valid 
10  August  1999, 
1200  UTC)  and 
MMS  metcm  (10 
August  1999, 
2200Z)  at  Profiler 
site. 
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Note  in  figure  18  is  the  +17°  K  virtual  temperature  difference  at  the 
surface.  The  fact  that  the  MMS  raob  is  from  mid-afternoon  explains  part 
of  the  difference.  More  specifically,  the  MMS  surface  temperature  would 
be  at  the  maximum  heating  portion  of  the  day  while  the  CAAM  metcm  is 
valid  when  there  is  maximum  cooling  at  the  surface.  Above  2  km,  the 
two  systems  agree  very  closely.  In  figure  19,  the  pressure  differences  are 
always  within  2  mb. 

3.5  Comparisons  of  a  CAAM  BFM  metcm  (from  a  Model 
Rim)  with  like  Data  Displaced  Either  Spatially  or 
Temporally  or  Both. 

Figures  20-23  illustrate  how  the  "stats"  program  can  be  utilized  to 
compare  the  u,  v,  Tv,  and  P  variables  between  two  CAAM  BFM  metcms. 
Here,  gnuplot  displays  data  from  CAAM  BFM  metcms  valid  at  Profiler 
site  and  Santa  Theresa,  NM  site  for  10  Aug  1999, 1200  UTC. 
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Figure  20.  U- 
component  deltas 
(m/ s)  between 
CAAM  metcms  at 
Profiler  site  (10 
August  1999, 1200 
UTC)  and  Santa 
Teresa,  NM  site  (10 
August  1999, 1200 
UTC). 


Figure  21.  V- 
component  deltas 
(m/s)  between 
CAAM  metcms  at 
Profiler  site  (10 
August  1999, 1200 
UTC)  and  Santa 
Teresa,  NM  (10 
August  1999, 1200 
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Figure  22.  Tv  deltas 
(deg  K)  between 
CAAM  metcms  at 
Profiler  site  (10 
August  1999, 1200 
UTC)  and  Santa 
Teresa,  NM  site  (10 
August  1999, 1200 
UTC). 
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Figure  23.  P  deltas 
(mb)  between  CAAM 
metcms  at  Profiler 
site  (10  August  1999, 
1200  UTC)  and  Santa 
Teresa,  NM  site  (10 
August  1999, 1200 


UTC). 
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As  seen  in  previous  cases,  as  one  gets  above  the  so-called  boundary  layer 
(at  about  1  km)  frictional  effects  become  less  important  and  stations  in 
relatively  close  proximity  in  time  and  space  show  fairly  good  agreement 
in  winds.  Temperature  gradients  will  decrease  with  height,  as  will  the 
effects  due  to  surface  type  differences.  In  addition,  the  local  geography 
becomes  less  and  less  of  a  factor.  These  factors  are  evident  here  as  well. 
In  figure  20,  above  2  km  the  u-component  differences  diminish 
significantly.  The  v-component  displays  the  same  characteristics. 
Temperature  deltas  are  virtually  zero  above  2  km.  The  pressure  deltas 
are  again  always  within  2  mb. 
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4.0  Conclusions 

The  stats  analysis  and  display  software  is  potentially  a  useful  tool  for 
determining  the  accuracy  of  profiler  data  when  used  as  input  to  the 
CAAM  BFM  mesoscale  model.  In  addition,  it  can  be  used  to  compare 
two  types  of  weather  balloon  data  formats,  the  MMS  and  those  such  as 
from  the  NWS.  Furthermore,  metcms  in  either  the  MMS  or  CAAM  BFM 
format  can  be  compared.  Thus  anywhere  the  CAAM  BFM  model  or 
MMS  balloon  data  is  being  analyzed  this  software  can  be  an  asset. 
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Acronyms 


AGL 

above  ground  level 

ARL 

U.S.  Army  Research  Laboratory 

CAAM  BFM 

Computer- Assisted  Artillery  Meteorology 
Battlescale  Forecast  Model 

metcm 

computer  meteorological  message 

MMS 

Meteorological  Measuring  Set 

MMS-P  (POC) 

Meteorological  Measuring  Set  -  Profiler 
(Proof-of-Concept) 

MSL 

mean  sea  level 

NWS 

National  Weather  Service 

raob 

radiosonde  observation 

rmse 

root  mean  square  error 

UTC 

Universal  Time  Coordinates 

Appendix  A 

Software  Module  Description 


i.  stats 

Description:  A  Bourne  shell  script.  The  user  types  "stats"  at  the  UNIX  prompt 
and  a  screen  will  be  displayed  giving  the  user  13  options: 

1)  Compare  MMS  raob,  CAAM  BFM  "out.bin"  data  and  display  results  with 
GNUPLOT 

2)  Compare  MMS  raob,  CAAM  BFM  "out.bin"  data  and  send  output  to  default 
printer 

3)  Compare  MMS  raob,  "NWS  format"  raob  and  display  results  with 
GNUPLOT 

4)  Compare  MMS  raob,  "NWS  format"  raob  and  send  output  to  default  printer 

5)  Compare  "NWS  format"  raob,  CAAM  BFM  "out.bin"  file  and  display  results 
with  GNUPLOT 

6)  Compare  "NWS  format"  raob,  CAAM  BFM  "out.bin"  file  and  send  output  to 
default  printer 

7)  Compare  2  "NWS  format"  raobs,  display  results  with  GNUPLOT 

8)  Compare  2  "NWS  format"  raobs,  send  results  to  default  printer 

9)  Compare  CAAM  BFM  metcm,  MMS  metcm  and  display  results  with 
GNUPLOT 

10)  Compare  CAAM  BFM  metcm,  MMS  metcm  and  send  results  to  default 
printer 

11)  Compare  2  CAAM  BFM  metcm's  and  display  results  with  GNUPLOT 

12)  Compare  2  CAAM  BFM  metcm's  and  send  results  to  default  printer 

13)  Exit  the  stats  program 

ii.  compare_raob_MMSraob.f 

Description:  Takes  "NWS  format"  raob,  MMS  format  raob,  performs  vertical 
interpolation  of  each  to  predefined  heights. 

Inputs: 

a.  NWS  format  raob 

b.  MMS  raob 

c.  "z_zm.file"  -  heights  at  which  the  pseudo-terrain-following  analysis  is 
produced. 

d.  "terrainl.in"  file  -  an  ascii  file  containing  gridded  elevation  data  in 
meters.  The  file  header  contains  the  center  point  latitude /longitude 
(deg),  the  number  of  grid  points  in  the  x-  and  y-directions,  and  the  grid 
spacing  (m). 


37 


Outputs: 

a.  "raob_compare_raobl" 

Contains:  latitude,  longitude,  height,  u-,  v-component,  ambient  and 
virtual  temperature,  and  pressure 

b.  "raob_compare_raob2" 

Contains  same  information 

iii.  compare_Taob_BFMout_bin.f 

Description:  For  comparison  of  a  "NWS  format"  raob  and  the  particular  forecast 
component  of  a  CAAM  BFM  out.bin  file.  Performs 

1)  Bilinear  interpolation  (i.e.  horizontal)  of  CAAM  BFM  gridded  data  to  get 
values  at  the  particular  latitude /longitude  of  the  raob  being  compared  with. 

2)  Vertical  interpolation  of  CAAM  BFM  data 

3)  Vertical  interpolation  of  NWS  format  raob 

Inputs: 

a.  "NWS  format"  raob 

b.  CAAM  BFM  "out.bin"  file 

c.  "z_zm.file" 

d.  "terrainl.in"  file 


Outputs: 

a.  "compare_raob" 

contains:  latitude,  longitude,  u-,  v-component,  ambient  and  virtual 
temperature,  and  pressure 

b.  "compare_model" 
contains:  same  information 

iv.  compare_raobl_raob2.f 

Description:  Takes  two  "NWS  format"  raobs  and  vertically  interpolates  each 
raob  to  predefined  heights.  In  order  to  calculate  virtual  temperature,  an 
intermediate  variable,  Q  (mixing  ratio),  is  defined  as  a  function  of  pressure  and 
dew  point  temperature. 

Inputs: 

a.  "z_zm.file" 

b.  "terrainl.in"  file 

c.  two  "NWS  format"  raobs 
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Outputs: 

a.  "raob_compare_raobl" 

contains:  latitude,  longitude,  u-,  v-component,  ambient  and  virtual 
temperature,  and  pressure 

b.  "raob_compare_raob2" 
contains:  same  information 

v.  parselntemetRaob.c 

Description:  By  searching  out  certain  keywords  and  positions  in  the  input  file, 
extracts  and  writes  to  the  output  file  the  variables:  Pressure,  height,  ambient  and 
dew  point  temperature,  wind  speed  and  direction. 

Input: 

a.  "NWS  format"  raob 
Output: 

a.  "Intemet_raob_reformatted" 

vi.  parseMMS.c 

Description:  Sifts  through  the  input  file  by  looking  for  keywords  and  certain 
positions  in  the  file  to  extract  the  variables:  pressure,  height,  ambient  and  dew 
point  temperature,  wind  speed  and  direction. 

Input:  MMS  raob 

Output:  "MMS_raob_reformatted" 

vii.  reform at_mms_me tcm . c 

Description:  Takes  a  MMS  computer  meteorological  message  which  has  multiple 
levels  of  data  per  line  of  wind  direction  and  speed,  virtual  temperature  and 
pressure. 

Input:  MMS  metcm 

Output:  "colformat_mms_metcm" 
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viii.  reformat_bfm_metcm_to_columns.c 

Description:  Takes  a  CAAM  BFM  computer  meteorological  message  which  has  3 
levels  per  line  of  wind  direction  and  speed,  virtual  temperature,  and  pressure. 

Input:  CAAM  BFM  metcm 

Output:  "colformat_bfm_metcm" 

Lx.  rd_raob_compare.c 

Description:  Assuming  "good"  values  for  the  variables,  u-  and  v-component, 
virtual  temperature,  and  pressure  at  matching  heights  in  the  raobs,  outputs  the 
differences  at  each  height  level. 

Inputs:  two  raobs  as  output  by  1  of  the  aforementioned  FORTRAN  programs 

Output:  "compare_raob_rmse"  containing  root  mean  square  errors  in  the  u-  and 
v-component,  virtual  temperature,  and  pressure 
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Appendix  B 
Program  Flow  Chart 


v.  stats  options  9  and  10 
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Appendix  C 

A  Sample  GNUPLOT  file 


In  this  particular  case  a  postscript  file  is  created.  To  change  it  to  screen  display  one  would 
remove  the  set  terminal  line  and  the  set  output  line. 


set  terminal  postscript 
#show  terminal 
set  autoscale 

set  output  "u_MMS_raob.ps" 
set  grid 

set  title  'Analysis' 

set  xlabel  'MMS  u-comp  -  raob  u-comp' 

#set  xtics  0,  2,  50 
#set  xrange[0  :  50] 

#set  xtics  -30,  2,  30 
#set  xrange[-30  :  30] 
set  ylabel  'Height  AGL  (m) '  2,0 
#set  ytics  0,  1000,  20000 
#set  yrange  [0  20000] 

#set  size  0.66,  0.8 

#set  label  1  "raob  Tv"  at  283.,  2250. 

#next  plot  diffs  in  u  and  v  respectively 

#plot  "mms_bfm_diffs"  using  3:2  notitle  with  linespoints  1  1 
#pause  -1 

plot  "compare_raob_delta_u_f iltered"  using  2:1  notitle  with 
linespoints  1  1 


NASA  MARSHALL  SPACE  FLT  CTR 
ATMOSPHERIC  SCIENCES  DIV 
E501 

ATTN  DR  FICHTL 
HUNTSVILLE  AL  35802 

NASA  SPACE  FLT  CTR 
ATMOSPHERIC  SCIENCES  DIV 
CODE  ED  41 1 
HUNTSVILLE  AL  35812 

US  ARMY  MISSILE  CMND 
AMSMIRDACAD 
ATTN  DR  PETERSON 
REDSTONE  ARSENAL  AL  35898-5242 

US  ARMY  MISSILE  CMND 
AMSMI RD  AS  SS 
ATTN  MR  H  F  ANDERSON 
REDSTONE  ARSENAL  AL  35898-5253 

US  ARMY  MISSILE  CMND 
AMSMI  RD  AS  SS 
ATTN  MR  B  WILLIAMS 
REDSTONE  ARSENAL  AL  35898-5253 

US  ARMY  MISSILE  CMND 
AMSMI  RDDESE 
ATTN  MR  GORDON  LILL  JR 
REDSTONE  ARSENAL  AL  35898-5245 

US  ARMY  MISSILE  CMND 
REDSTONE  SCI  INFO  CTR 
AMSMI  RDCSR  DOC 
REDSTONE  ARSENAL  AL  35898-5241 

US  ARMY  MISSILE  CMND 
AMSMI 

REDSTONE  ARSENAL  AL  35898-5253 

PACIFIC  MISSILE  TEST  CTR 
GEOPHYSICS  DIV 
ATTN  CODE  3250 
POINT  MUGU  CA  93042-5000 

ATMOSPHERIC  PROPAGATION  BRANCH 
SPAWARSYSCEN  SAN  DIEGO  D858 
49170  PROPAGATION  PATH 
SAN  DIEGO  CA  92152-7385 

METEOROLOGIST  IN  CHARGE 
KWAJALEIN  MISSILE  RANGE 
PO  BOX  67 

APO  SAN  FRANCISCO  CA  96555 


NCAR  LIBRARY  SERIALS 
NATL  CTR  FOR  ATMOS  RSCH 
PO  BOX  3000 
BOULDER  CO  80307-3000 

HEADQUARTERS  DEPT  OF  ARMY 
DAMIPOI 
ATTN  LEE  PAGE 
WASHINGTON  DC  20310-1067 

MIL  ASST  FOR  ENV  SCI  OFC 
OF  THE  UNDERSEC  OF  DEFNS 
FOR  RSCH  &  ENGR  R&AT  E  LS 
PENTAGON  ROOM  3D129 
WASHINGTON  DC  20301-3080 

DEANRMD 
ATTN  DR  GOMEZ 
WASHINGTON  DC  20314 

US  ARMY  INFANTRY 
ATSHCDCSOR 
ATTN  DR  E  DUTOIT 
FT  BENNING  GA  30905-5090 

HQ  AFWA/DNX 

106  PEACEKEEPER  DR  STE  2N3 

OFFUTT  AFB  NE  68113-4039 

PHILLIPS  LABORATORY 
PLLYP 

ATTN  MR  CHISHOLM 
HANSCOM  AFB  MA  01731-5000 

PHILLIPS  LABORATORY 
PL  LYP  3 

HANSCOM  AFB  MA  01731-5000 

AFRL/VSBL 
29  RANDOLPH  RD 
HANSCOM  AFB  MA  01731 

US  ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY 

ATTN  MR  H  COHEN 
APG  MD  21005-5071 

US  ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY  AT 

ATTN  MR  CAMPBELL 
APG  MD  21005-5071 

US  ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY  CR 
ATTN  MR  MARCHET 
APG  MD  21005-5071 


ARL  CHEMICAL  BIOLOGY 
NUC  EFFECTS  DIV 
AMSRL  SL  CO 
APG  MD  21010-5423 

US  ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMSXY 

APG  MD  21005-5071 

ARMY  RESEARCH  LABORATORY 
AMSRL  D 

2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  OP  Cl  SD  TL 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  LL 
ADELPHI  MD  20703-1197 

ARMY  RESEARCH  LABORATORY 
AMSRL  SSSH 
ATTN  DR  SZTANKAY 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl 
ATTN  J  GANTT 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1197 

ARMY  RESEARCH  LABORATORY 
AMSRL 

2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1145 

NATIONAL  SECURITY  AGCY  W21 
ATTN  DR  LONGBOTHUM 
9800  SAVAGE  ROAD 
FT  GEORGE  G  MEADE  MD  20755-6000 

US  ARMY  RSRC  OFC 
ATTN  AMXRO  GS  DR  BACH 
PO  BOX  12211 
RTP  NC  27009 

DR  JERRY  DAVIS 
NCSU 

PO  BOX  8208 
RALEIGH  NC  27650-8208 


US  ARMY  CECRL 
CECRLGP 
ATTN  DR  DETSCH 
HANOVER  NH  03755-1290 

USARMYARDEC 
SMCAR IMI I  BLDG  59 
DOVER  NJ  07806-5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MT  DA  L  3 
DUGWAY  UT  84022-5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MTM 
ATTN  MR  BOWERS 
DUGWAY  UT  84022-5000 

DEPT  OF  THE  AIR  FORCE 

OL  A  2D  WEATHER  SQUAD  MAC 

HOLLOMAN  AFB  NM  88330-5000 

PL  WE 

KIRTLAND  AFB  NM  87118-6008 

USAF  ROME  LAB  TECH 
CORRIDOR  W  STE  262  RL  SUL 
26  ELECTR  PKWY  BLD 106 
GRIFFISS  AFB  NY  13441-4514 

AFMCDOW 

WRIGHT  PATTERSON  AFB  OH  45433-5000 

US  ARMY  HELD  ARTILLERY  SCHOOL 

ATSFTSMTA 

FT  SILL  OK  73503-5600 

US  ARMY  FOREIGN  SCI  TECH  CTR 
CM 

220  7TH  STREET  NE 
CHARLOTTESVILLE  VA  22448-5000 

NAVAL  SURFACE  WEAPONS  CTR 
CODE  G63 

DAHLGREN  VA  22448-5000 

US  ARMY  OEC 

CSTE  EFS 

PARK  CENTER  IV 

4501  FORD  AVE 

ALEXANDRIA  VA  22302-1458 

US  ARMY  CORPS  OF  ENGRS 
ENGR  TOPOGRAPHICS  LAB 
ETLGS  LB 

FT  BELVOIRVA  22060 
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US  ARMY  TOPO  ENGR  CTR 
CETEC  ZC  1 

FT  BELVOIR  VA  22060-5546 

SCI  AND  TECHNOLOGY 
101  RESEARCH  DRIVE 
HAMPTON  VA  23666-1340 

US  ARMY  NUCLEAR  CML  AGCY 
MONA  ZB  BLDG  2073 
SPRINGFIELD  VA  22150-3198 

USATRADOC 

ATCDFA 

FT  MONROE  VA  23651-5170 

ATRCWSSR 
WSMR  NM  88002-5502 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  S 
COMP  &  INFO  SCI  DIR 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  E 
COMP  &  INFO  SCI  DIR 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  W 
COMP  &  INFO  SCI  DIR 
WSMR  NM  88002-5501 

DTIC 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BELVOIR  VA  22060-6218 

US  ARMY  MISSILE  CMND 
AMSMI 

REDSTONE  ARSENAL  AL  35898-5243 

US  ARMY  DUGWAY  PROVING  GRD 
STEDP3 

DUGWAY  UT  84022-5000 

USTRADOC 

ATCDFA 

FT  MONROE  VA  23651-5170 

WSMR  TECH  LIBRARY  BR 
STEWS  IM  IT 
WSMR  NM  88002 

US  ARMY  RESEARCH  LAB 
AMSRL  D  DR  D  SMITH 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1197 


USARMYCECOM 
INFORMATION  &  INTELLIGENCE 
WARFARE  DIRECTORATE 
ATTN  AMSEL  RD IW  IP 
FORT  MONMOUTH  NJ  07703-5211 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  EA 
S  KIRBY 

COMP  &  INFO  SCI  DIR 
WSMR  NM  88002-5501 


Record  copy 
Total 
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